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Abstract 



The paper is concerned with mechanical systems which are controlled by implementing 
a number of time-dependent, frictionless holonomic constraints. The main novelty is due 
to the presence of additional non-holonomic constraints. We develop a general framework 
£>. ■ to analyze these problems, deriving the equations of motion and studying the continuity 

I/"") , properties of the "control-to-trajectory" maps. Various geometric characterizations are 

provided, in order that the equations be affine w.r.t. the time derivative of the control. In 
this case the system is fit for jumps, and the evolution is well defined also in connection 
with discontinuous control functions. The classical Roller Racer provides an example where 
the non-afhne dependence of the equations on the derivative of the control is due only to 
the non-holonomic constraint. This is a case where the presence of quadratic terms in the 
equations can be used for controllability purposes. 

1 Introduction 



The control of mechanical systems provides a rich area for mathematical investigation. In a 
commonly adopted model [2, 3, 12, 22], the controller can modify the time evolution of the 
system by applying additional external forces. This leads to a control problem in standard 
form, where vector fields governing the state variables depend continuously on the control 
function. 

In an alternative model, also physically meaningful, the controller acts on the system by 
directly assigning the values of some of the coordinates, as functions of time. Stated in a 
more intrinsic fashion, this means that the controller assigns the leaves of a foliation of the 
configuration manifold as functions of time. Here the basic framework consists of a manifold 
with coordinates q = (q 1 , . . . , q N ,q N+1 , ■ ■ ■ , q N+M ). We assume that the values of the last M 
coordinates can be prescribed: q N+1 = u 1 ^), q N+M = u M (t). Instead of forces, these 
controls thus take the form of time-dependent holonomic constraints. The evolution of the 
remaining free coordinates q 1 , . . . ,q N is then 
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determined by a control system where the right hand side depends not only on the control 
itself, but also (linearly or quadratically) on the time derivative of the control function. 

This alternative point of view was introduced, independently, in [9] and in [19]. A con- 
siderable amount of literature is now available on mechanical systems controlled by active 
constraints. The form of the basic equations of motion, in relation with geometrical proper- 
ties of the system, was studied in [1, 9, 13, 17, 23, 24]. Differential equations of impulsive 
nature, where the right hand side contains a measure, such as the distributional derivative 
of a discontinuous control, were considered in [5, 6, 7, 18, 20, 25, 28]. These works focus, 
in particular, on the continuity of the "control-to-trajectory" map, in various topologies. In 
addition, problems of stabilization and of optimal control for this kind of "hyper-impulsive" 
mechanical systems were studied in [1, 8, 16] and in [6, 7, 10, 11], respectively. See also [4] for 
a survey. 

The goal of the present paper is to extend this theory to the case where the system 
is subject to some additional non-holonomic constraints. In this case, the velocity vector 
q = (q 1 , . . . , q N + M ~j satisfies an additional set of v linear relations 

N+M 

]T Jtd = o fc = i,. ..,!/. 

1=1 

A major focus of our analysis is on the form of the resulting equations, In general, the 
right hand side of the evolution equations turns out to be a quadratic polynomial of the time 
derivative u( ). Therefore, the evolution problem is well posed as soon as the control function 
satisfies u(-) <G W 1,2 . There are, however, important cases where the right hand side is an 
affine function of u(-). Our main result in this direction, Theorem 6.1, yields a number of 
equivalent analytic and geometric conditions for this to happen. In the positive case, the 
"control-to-trajectory" map can be extended by continuity to a larger family of (possibly 
discontinuous) control functions. Following [9], we then say that the system is "fit for jumps". 
Section 7 provides an additional geometric characterization of this property. Let 

A u = {(q\...,q N+M ); q N+a = u a , 1 < a < M} (1.1) 

be the leaf of the foliation corresponding to the control value u. Roughly speaking, we show 
that the system is "fit for jumps" if and only if the "infinitesimal non-holonomic distance" 
between leaves of the foliation remains constant, along directions compatible with the non- 
holonomic constraint. 

For sake of illustration, in the last section we show how the present framework applies to 
the control of the Roller Racer. An interesting feature of this example is that, without the 
non-holonomic constraint, the equations of motion would be affine w.r.t. the time derivative it 
of the control function. However, the additional non-holonomic constraint renders the system 
not fit for jumps. It is indeed the presence of a quadratic term in the derivative of the control 
that makes forward motion possible. 

2 Non holonomic systems with active constraints as controls 

Let N, M, v be positive integers such that v < N + M. Let Q be an (N + M)-dimensional 
differential manifold, which will be regarded as the space of configurations of a mechanical 
system. 

Let r be a distribution on Q, i.e. a vector sub-bundle of the tangent space TQ. Throughout 
the following, we consider trajectories t i-> q(i) G Q of the mechanical system which are 
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continuously differentiable and satisfy the geometric constraint 

q(t) G r q(t) . (2.1) 

We do not assume T to be integrable, so that in general (2.1) is a non holonomic constraint. 

In our model, the system will be controlled by means of an active (holonomic, time- 
dependent) constraint. To describe this constraint, let an M-dimensional differential manifold 
U be given, together with a submersion 

ir : Q\->IA. (2.2) 

The fibers 7r _1 (u) C Q will be regarded as the states of the active constraint. The set of all 
these fibers can be identified with the control manifold U. 

Let X be a time interval and let u : X h-> U be a continuously differentiable map. We say 
that a trajectory q : X — > Q agrees with the control u(-) if 

vroq(t) = u (t) for all teX. (2.3) 

For each q G Q, consider the subspace of the tangent space at q given by 

A q = ker T q 7r . 

Here T q 7r denotes the linear tangent map between the tangent spaces T q Q and T n ^U. Clearly, 
A is the (holonomic) distribution whose integral manifolds are precisely the fibers 7r -1 (u). 



2.1 General setting 

1) The manifold Q is endowed with a Riemannian metric g = g q [-, •], the so-called kinetic 
metric, which defines the kinetic energy T. More precisely, for each q G Q and v G T q Q 
one has 

T(q,v) = ^g q [v,v]. (2.4) 

We shall use the notation v i-> g q (v) to denote the isomorphism from T q Q to T*Q 
induced by the scalar product g q [-, •]. Namely, for every v G T q Q, the 1-form g q (v) is 
defined by 

(gq(v),w) = gq[v,w] for all w G T q Q, (2.5) 

where (■, •) is the natural duality between the tangent space T q Q and the cotangent 
space T* Q. 

If q G Q and W C T^, W 1 - denotes the subspace of T q consisting of all vectors that are 
orthogonal to every vector in W: 

W 1 - = {v G T q I g q [v, w] = for all w G W} . 

For a given distribution E C TQ, the orthogonal distribution E 1 - C TQ is defined by 
setting E^r = (T^)" 1 , for every q G Q. 

2) Throughout the following, we shall assume that the holonomic distribution A and the 
non-holonomic distribution T satisfy the tranversality condition 

Aq + r q = T q Q for all q G Q . (2.6) 

Notice that this is equivalent to 

A 1 nr 1 = {o}, (2.7) 

and implies v < N. 
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3) The mechanical system is subject to forces. In the Hamiltonian formalism, these are 
represented by vertical vector fields on the cotangent bundle T*Q. We recall that, 
in a natural system of coordinates (q,p), the fact that F is vertical means that its q- 
component is zero, namely F = J2i=i M 

4) The constraints (2.1) and (2.3) are dynamically implemented by reaction forces obeying 

D'Alembert condition: If t ^ q(£) is a trajectory which satisfies both (2.1) and 
(2.3), and R(i) is the constraint reaction at a time t, then 

R(t) G ker(A q(t) nr q(t) ). (2.8) 

In other words, regarding the reaction force R(t) as an element of the cotangent space 
T* (t) Q, one has 

(R(t), v) = for all v G A q(t) nr q(() C T q(t) Q. (2.9) 



2.2 Equations of motion 

For each q € Q, we shall use g" 1 to denote the inverse of the isomorphism g q at (2.5) 
Moreover we define the scalar product on the cotangent space g q 1 [-, ■ 
by setting 



T q Q x T^Q i y M 



gq^.P] = gqfeq^P), gq^P)] 



for all (p, p) 



€ T*QxT*Q. 



For every q G Q, we shall use "H(q, •) to denote the Legendre transform of the map v 
T(q, v), so that 



«(q,P) = W[P,P] (=T(q,g- 1 (p))) for all p € T q Q . 



(2.10) 



The map "H : T*Q — > iR will be called the Hamiltonian corresponding to the kinetic energy T ■ 
Let (g) and (n) be local coordinates on Q and ZY, respectively, such that the domain of 
the chart (q) is mapped by it into the domain of the chart (u). Let (q,p) the natural local 
coordinates on T*Q corresponding to the coordinates (q). From Nonholonomic Mechanics it 
follows that, given a smooth control t i-> u(t) (here regarded as a time-dependent holonomic 
constraint), the corresponding motion t i-» (q,p)(t) on T*Q verifies the relations 



p(f) 

7T O g(t) 



dp 



(Q(t),p(t)), 



dH 

— (q{t),p(t)) + F(t,q(t),p(t)) + i?(t) , 



«(t), 

p(t) e ff g (t)(r g ( t ))> 

-R(t) G kerA 9(t) +kerr g(t) 



(2.11) 



We have used H, F, R, and g to denote the local expressions of H, F, R, and g, respectively. 
For simplicity, we use the same notation for the distributions V, A on Q and their local 
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expressions in coordinates. Moreover, we use the notation v i-> g q (v) to denote the local 
expression of the isomorphism v i— >■ g q (v). 

The first two equations in (2.11) are the dynamical equations written in Hamiltonian form. 
The first one yields the inverse of the Legendre transform. The third equation represents 
the (holonomic) control-constraint. Relying on the fact that it is a submersion, we shall 
always choose local coordinates (q r ) and (u Q ) such that q N+a = u a , for all a = 1, . . . , M. 
We thus regard this equation as prescribing a priori the evolution of the last M coordinates: 
q N+1 , . . . , q N+M . The fourth relation in (2.11) is the Hamiltonian version of the non holonomic 
constraint (2.1). The fifth relation is clearly equivalent to (2.8), i.e. it represents d'Alembert's 
condition. 

Remark 2.1 Although a global, intrinsic formulation of these equations can be given (see 
[19]) we are here mainly interested in their local coordinate- wise expression. Indeed, a major 
goal of our analysis is to understand the functional dependence on the time derivative u of 
the control. 



The special case where no active constraints are present can be obtained by taking A = T*Q, 
i.e. ker(A) = {0}. In this case, (2.11) reduces to the standard Hamiltonian version of the 
dynamical equations with non-holonomic constraints, namely 

dH 



pit) 



dH 
dq 



(q(t),p(t))+F(t,q(t),p(t))+R(t), 



(2.12) 



p(t) e 0,(t)(r,(t)), 

{ R(t) G kerr, w , 



3 Orthogonal decompositions of the tangent and the cotangent 
bundles 

To derive a set of equations describing the constrained motion, it will be convenient to de- 
compose both the tangent bundle TQ and the cotangent bundle T*Q as direct sums of three 
suitable vector sub-bundles. We recall that Q is a manifold of dimension N + M, while A and 
T are distributions on Q, having dimensions N and N + M — u, respectively. (In view of the 
transversality condition (2.6) one has v < N.) 

3.1 Tangent bundle 

Definition 3.1 For every q G Q, we define the following three subspaces o/T q Q. 

(T q Q) 7 = A q nr q , (T q Q) 77 = r£, (r q Q) m = (A, n r,)- 1 n r q . (3.1) 

Proposition 3.1 For each q G Q, the three subspaces in (3.1) are mutually orthogonal and 
span the entire tangent space, namely 

TqQ = (T q Q)/ © (TqQ) n © (T q Q)/// , (3.2) 
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Anr=(TQ) I 



rn(AnrV-=cTQ) 




Figure 1: The orthogonal decomposition of TQ. 

Moreover, 

{TqQ)l © {TqQ)lII = T. (3.3) 

If the transversality condition (2.6) holds, then the above subspaces have dimensions 

dim(T q Q) 7 = N-v, dim(T q Q) 77 = v, dim(T q Q) m = M. (3.4) 



Proof. The orthogonality of the subspaces (T q Q)j and (T q Q)/j is immediately clear form 
the definitions. Observing that 

(jTqQ)j © (TqQ)jj^ ± = ((T q Q) / ) ± n((T q Q) // ) ± = (AqHr/flTq = (T q Q) m , 
we obtain the orthogonal decomposition (3.2). In particular, this implies 

{TqQ)i © {TqQ)in = ((TqQ)//) 1 " = r q . 

Finally, if (2.6) holds, then 

dim(A q nr q ) = dim(A q ) + dim(r q )-dim(T q Q) = N + (N + M - v) - (N + M) = N-v. 

Moreover, dim((r q Q)//) = dim(T q Q) - dim(T q ) = v. The last equality in (3.4) now follows 
from (3.2). □ 



by 



For J G {I, II, III}, the orthogonal projection onto the subspace {TQ)j will be denoted 



Vj : T q Q h+ (T q Q)j 



(3.5) 



Remark 3.1 Under the assumption (2.6), the third subspace in the decomposition (3.2) can 
be equivalently written as 

(T q Q)m = V r (A^), 



where V T denotes the orthogonal projection on the space V. Indeed, for any vector v G T q <2 
one has v G V T {^) iff 

v g r q n(A^ + r^) = r q n(A q nr q )^. 

3.2 Cotangent bundle 

Thanks to the isomorphism g : TQ T*Q defined at (2.5), one can use (3.2) to obtain a 
similar decomposition of the cotangent bundle as the direct sum of three vector bundles: 

T*Q = (T*Q)/ © (T*Q) n © (T*Q) in , (3.6) 

where, for J G {I, II, III}, 

(T^Q)j = ^((TqQJj). (3.7) 

We denote by Vj : T q Q — > (T q Q)j the orthogonal projection w.r.t. the metric g _1 . The above 
construction yields 

V} = goVjog- 1 . (3.8) 

Remark 3.2 By (3.8), the orthogonal decompositions of the tangent and cotangent bundles 
at (3.2) and (3.6) have the following property. Recalling (2.5), assume that p = g(q), so that 
q = g-^p). Then, for J G {/, II, III}, 

Pj = v*M = g(r.M)), qj = VM) = g-\v*j( P )). (3.9) 

In other words, the J-th component of p depends only on the J-th component of q, and 
viceversa. 

In view of Proposition 3.1 one obtains: 
Proposition 3.2 If the transversality condition (2.6) holds, then 

dim((T*Q) / ) = N — v , dim((T*Q) JJ ) = u, dim((T*Q) m ) = M. (3.10) 

Moreover, for every q G Q the three subspaces (T*Q)i, (T*Q)u, (T*Q)iu are pairwise orthog- 
onal (w.r.t. the metric g^ 1 ). In particular, (T*Q)u = ker(r q ). 

Next, we observe that by (2.6) the differential of the submersion ir : Q h-> U is one-to-one 
when restricted to (T q Q)jjj. More precisely, Dtt : (TqQ)jjj >->■ T n ^U is a bijective linear 
map. Its inverse will be denoted as 

h:T ff(q) W ^ (T q Q) m . (3.11) 

For any v G T W ^U, by h(v) we thus denote the unique vector in (T q Q)/jj such that 

L>7r-(h(v)) = v. (3.12) 

Moreover, we write 

k(v) = g q (h(v)) G (T*Q) in . (3.13) 
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We claim that the vector h(v) defined above can be characterized as the unique vector w G T q 
where the following constrained minimum is attained: 

h(v) = argmin g q [z,z]. (3.14) 

ZGTq , Dn Z,=V 

Indeed, since the Riemann metric g q is positive definite, the right hand side of (3.14) is well 
defined. Consider any vector w G T q Q such that 

w G T q , Dtt- w = v. (3.15) 

The above equalities imply 

w = h(v) + w/ 

for some wj G (T q Q)j. Since the vectors h(v) G (TqQ)/// and wj are orthogonal, we have 

gq[w,w] = gq[h( V ), h(v)] +gq[w/,W/] > g q [h( V ), h(v)]. 

This proves our claim. 

Remark 3.3 Let ueW and choose q G Q such that u = vr(q). Identifying the tangent spaces 
T U U w T q Q/A q , the inner product 

(v,v')u = gq[h(v), h(v')] 

can be seen as a Riemann metric on the space of leaves. In general, this metric is not canon- 
ically defined, because it depends on the choice of a particular point q on each leaf. From a 
more precise result proved in Section 7 it will follow that, if this metric is independent of the 
choice of q, then the system is "fit for jumps" . Namely, the right hand sides of the dynamical 
equations (4.8) are affine functions of the time derivative u of the control. As a consequence, 
trajectories of the system can be meaningfully defined also in connection with discontinuous 
control functions [5, 6]. 



4 A closed system of control equations 

In the original formulation (2.11), the motion is characterized in terms of a family of ODEs 
coupled with a set of constraints. Relying on the decompositions at (3.2) and (3.6), we will 
show that the motion can be described by a system of differential equations. 

The following result, providing different ways to express the constraint (2.1), is straight- 
forward. 

Lemma 4.1 Let t i— >■ q(i) G Q be a C 1 map. Using the same notation as in (3.9), the non- 
holonomic constraint q(i) G r q ( t ) can be expressed in any of the following equivalent forms: 

p(i)eg(r q(t )) q//W = o P //(t) = o. (4.1) 

On a given, natural chart (q,p) on T*Q, and for J G {/, //, III}, we denote by Pj the (in- 
dependent) matrix representing the projection Vj. Let g = (g r , s ) be the matrix representing 
the Ricmannian metric g in the (/-coordinates. In turn, the inverse matrix g^ 1 = (g r,s ) 
represents the metric g^ 1 on the cotangent space. In this coordinate system, the components 
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of the force and of the constraint reaction will be denoted as Fj = Pj F and Rj = PJ R. 
The coordinate representations of the linear maps h and k introduced at (3.11)-(3.13) will be 
denoted by h and k, respectively. 

Throughout the following, we use a local system of adapted coordinates, so that q N+a = u a 
for a = 1, . . . , M. If the holonomic active constraint is satisfied, by (3.9) and (3.12)-(3.13) 
this implies the identity 



PnMt)) = h(u(t)), 



PiiMt)) = Hu(t)). 



(4.2) 



In terms of the orthogonal decompositions (3.2) and (3.6), the non-holonomic and the 
active holonomic constraints yield the system of equations 



qil(t) = 
qm(t)= h{u{t)) 



( 



which holds if and only if 



V 



Pii(t) = o 
pm(t)= k(u(t)) 



(4.3) 



To complete the description of the motion, it remains to derive the equations for qj and 
pi. By (3.9) it follows 



qi = Pj{q) = Pi{g-\ P )) = g-\P!{ P )). 
Moreover, differentiating pj = Pj(p) w.r.t. time and using (2.11), we obtain 



(4.4) 



Pi 



dP? 



Oq 
~8q~ 



L -q){p)+PKp) 



(g-\pi) + h(uj) ) (pj + pni) + PI (-^(q(t),p(t)) + F(t, q(t),p(t)) + R(t) 



Oq 



^ • g-\pi + k(u)fj (pj + k(u)) - \p*j (%^-Ipi + Ku) , Pi + k(u)]j + Fj . 

(4-5) 

Indeed, the fifth relation in (2.11), i.e. d'Alembert's condition, implies Pf(R(t)) = 0. 

The first two terms on the right hand side of (4.5) can be written in a simpler form, 
introducing the bilinear ((/-dependent, possibly not symmetric), ]R N+M - valued map 



(p,p) ^ 0i\p,p] = 



dp; 

dq 



1 



g-\p)){p)--Pl 



' dg' 1 
dq 



\P,P})- 



Setting 

T[pj,«] = e/[pj,fc(«)] + ej[fc(«), pi], V[u,u] = 

the equations (4.4) and (4.5) can be written in the form 1 

qi(t) = g~ l {pi), 

< 

( Pi(t) = 0i]pi,Pi] + *[«,«] 

We can now state the main result of this section. 



(4.6) 



(4.7) 



(4.8) 



x We recall that gj(t) = P/(4f ) while pi(t) = 4(P/(p)) 
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Theorem 4.1 Let 1 be a time interval and let u : I — > M be a C control function. Let 
q : X i— >■ IR N be a C 1 path and set p{t) = g(q(t)). Then the path (q,p) :T^T*Q satisfies the 
non holonomic equations (2.11) (with constraints as controls) if and only if the following two 
conditions are satisfied. 

(i) For J € {I, II, III}, the components qj = Pj(q(t)) and p,j{t) = Pj(p(t)) satisfy the 
equations in (4-3) and (4-8). 

(ii) At some time to E Z one has q N+a (to) = u a {to) for all a = 1, ... , M. 



Proof. 1. By the previous analysis, if all relations in (2.11) are satisfied, then the equations 
(4.3) and (4.8) hold. Moreover, the condition it o q(t) = u(t) implies (ii), for all times 

2. Next, assume that all the equations in (4.3) and (4.8) hold, and 7r(g(to)) = u(to) at 
some time to £ By (4-1), the equation qu = implies the fourth relation in (2.11). 
For any t G I, if vr(g(t)) — u(t) = then 

J t (*(q(t)) -«(*)) = Dit ■ q(t) - ii(t) = Dtt ■ q in {t) - u(t) = 0. 

From the assumption w(q(to)) — u(to) = and the regularity of the coefficients of the equation, 
we conclude that ir(q(t)) — u(t) = for all t El. Hence the third relation in (2.11) holds. 
The equations for q in (4.3) and (4.8) imply 

OH 

q = qi + qii + qin = g^ipi) + h(u) = ff _1 (pj + pn + pm) = -^—(q,p), 



dp 



This yields the first equation in (2.11). 
Finally, by defining 



f) J-f 

R(t) = P+-^(q,p)-F(t,q(t),p(t)), (4.9) 

the second equation in (2.11) is clearly satisfied. It remains to check the fifth equation in (2.11), 
namely R(t) G ker A g (t)+kerT q (t). Since by construction ker A +kerT = (T*Q) n + (T*Q) in , 
this is equivalent to proving that Rj = 0. Using the second equation in (4.8) we find 

Rj = Pf{R) = P*(p) + P*^L(q,p^- Fl 
= p I -^L.q) { p ) + P I {^(q, P ))-F I 

= pi- • g-\pi + m)) (pi + m) + \pi (%V + > pi + k ^)}) - Fi 



= 0. 



□ 
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5 The equations in A-adapted coordinates 

Let q = (g 1 , . . . , q N+M ) be A-adapted coordinates on a open set O C Q , so that 

r 9 5 i 



The main goal of this section is to write the equations of motion explicitly in terms of these 
state coordinates, together with adjoint coordinates £j, 1 < j < N + M corresponding to 
suitable bases of the cotangent bundle T*Q, decomposed as in (3.6). It will turn out that the 
relevant equations will involve only the 2N — v variables q 1 , . . . , q N , £i, . . . , £n-i>- 

Consider a family {Vi, . . . , ^Vn+m} of smooth, linearly independent vector fields on Q, 
such that, at least on the open set O, 

(T q Q) 7 = apan{Vi(q),...,V JV _ I ,(q)}, 

(T q Q)/j = apan{v JV _ I , + i(q),...,V; V (q)}, (5.1) 

(r q Q)/7/ = span|Vjv+i(q),... ,V^ +M (q)}- 

Throughout the following, we assume that the vectors {Vi, ... , Vjv-i/}, which generate 
(TqQ)i, are mutually orthogonal, i.e. 

g[V r ,V s ] =0 for all r, s € {1, . . . ,N — is}, r^s. (5.2) 

In addition, for % = 1, . . . , N + M, we define the basis of cotangent vectors 

* " TO (5 ' 3) 

By (3.7), this yields 

(T*Q)i = S pan{fii(q),...,fljv-,(q)}, 
(r*Q)// = span (q),...,^(q)}, 

(r*Q)m = span{n N+1 (<i), . . . ,Q, N+M (q)] ■ 

By (5.2) and (5.3), the differential forms {fii, ... , Qn-v} are mutually orthogonal with re- 
spect to the metric g _1 , namely 

g _1 [fi r , n„] =0 for all r,s E {1, . . . ,N — u}, r^s. (5.4) 
Moreover, the basis jfii(q), . . . , fi;v-i/(q)} is dual to the basis jVi(q), . . . , Vjv-j/(q)}, i.e. 

ft r (q) (V s (q)) = 5 riS for all r, s = 1, . . . , JV - i/. 

This choice of orthogonal bases makes it easy to compute the projections Vi and V\. Indeed, 
for any tangent vector w and any cotangent vector p one has 



n 



N-v 



g _1 [p,n 



^?(p) = e ( 5 - 6 ) 

In the following, to simplify notation, whenever repeated indices taking values from 1 to 
N + M are summed, the summation symbol will be omitted. On the other hand, summations 
ranging over a smaller set of indices will be explicitly written. As before, we let g = (g r , s ) be 
the matrix representing the Ricmannian metric g in the ^-coordinates. In turn, the inverse 
matrix g^ 1 = (g r,s ) represents the metric g _1 on the cotangent space. 

For t = 1, ... ,N + M, let V t x , ... , V e N+M be the g-components of V e , so that V t = V l i -J^.. 
If the (column) vector w = (w s ) € M N+M yields the coordinate representation of w, then by 

(5.5) the projected vector Pi(w) is a column vector with coordinates (Piw) s = (Pi)f.w r , 
where the (N + M) x (N + M) matrix Pj is defined by 

(Pj )' = y Sr ' fc L 1 r,s = l,...,N + M . 

Similarly, let fi rj i, . . . ,£l r ,N+M be the components of fi r , so that fi r = £l rtS dq s . If the 
(row) vector p G M M+N yields the coordinate representation of the covector p, then by 

(5.6) the projected the vector Pj(p) has coordinates given by (pPj) s = p r (Pj) r s , where the 
(N + M) x (N + M) matrix Pf is defined by 

( p*Y ^ 9 r,k ^e,k Me,s (c . 7] 

i l)s " h ' ( } 

In order to write the second set of equations in (4.3), we need an explicit expression of the 
((/-dependent) matrices h and k. Let us define the (N + M) x M matrix Vni and the M x M 
matrix Vm 111 by setting 



Vni = (V N+Ct r ) , V in in = (v N+a N+ ?) 



Here and in the sequel, Greek indices such as a, (3 range from 1 to M, while Latin indices such 
as r, s range from 1 to N + M. Notice that, by (5.1), the columns of Vm span (T„Q)jjj. 

Recalling the definitions (3.12)-(3.13), and the identities q N+a = u a , it is easy to check 
that the injective linear map h : TU W ^ i-> (T CL Q) III C T q Q is represented by the (N + M) x M 
matrix 

h = V IH ■ (Vm 111 )' 1 ■ 
In turn, the linear map k(-) = g(h(-)) is represented by the (N + M) x M matrix 



k = g ■ h = g-Vm- (Vm 111 ) 



-l 



Using this particular system of coordinates, the equations of motion (4.3), (4.8) for the 
non-holonomic system with active constraints can be written in the following form. 

Proposition 5.1 Let (q,p,u)(-) be as in Theorem 4-1- Moreover, let (£, 77, A) be the compo- 
nents of p = pidq 1 w.r.t. the frame {fii, . . . , Qn+m}, so that 

N-v N N+M 

P(t) = E E Vi(t)Cli+ E MWe- 

t=l l=N-v+l i=N+l 
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Then the curve t 1-> (q(t) , p(t) , u(t)) provides a solution to the system (2.11) if and only if the 
curve t i-> (q(t), £(t), r](t), X(t), u(t)) is a solution of 



V/ 



M 
a=l 



q N+a = u a 



r=l,...,N, 

a = 1,...,M , 
m = 1, . . . , N — v , 



(5.8) 



im = e m [s, Z} + f m [Z,u\ + V m [u,u\ + F m 

v = o, 

A = {Vnif-k-u, 

where the superscript denotes transposition and, recalling (4-6)-(4- 7), for every m = 1, . . . , N— 
v we set 



N 



= Ea,b=i otfUb, with e a J = ( [ej]p a n r> «J2 8)6 - £ -^V' s ^, 6 ) v m \ 



N 



dq N+a Q q j 



^ m [w,w] = y a> l3,mW a wP, With y a ,i3,m = ^a,/3,iV m l , 



F — FV 1 



(5.9) 



Indeed, (5.9) provides the representation of the maps 0/,T, and ^ in (4.6)-(4.7), with the 
present choice of coordinates. 



6 Systems which are "fit for jumps" 

We now examine the connections between the following properties: 

i) The continuity of the control-to-trajectory map u(-) i— > (g(-), £>(•)). 

ii) The vanishing of the "centrifugal" term u] in (4.7). 

hi) The invariance of the distribution r n (r D A)- 1 in T-constrained inertial motions. 

We recall that, if it] = 0, then the equations of motion (5.8) are affine w.r.t. the 
time derivative ii of the control function. In this case, the control-to-trajectory map u(-) i— > 
(q(-),p(-)) can be continuously extended to a family of possibly discontinuous control functions 
u(-) with bounded variation [5]. As in [9], we then say that the system is fit for jumps. On the 
other hand, if the quadratic term 'J/f-ii,^] does not vanish, this means that small vibrations of 
the control produce a centrifugal force. As shown in [8], this feature can be exploited as an 
additional way to control or stabilize the system. 

As in Section 2, T will denote the kinetic energy associated with the metric g, defined at 
(2.4), while T denotes a non-holonomic distribution on the manifold Q. 
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Definition 6.1 Let X be a time interval. A C 2 map q : X — >■ Q will be called a free inertial 
motion if, in any set of coordinates, its local expression q(-) provides a solution to the Euler- 
Lagrange equations 

dtm~^ = °- ((U) 

Definition 6.2 A C 2 map q : X — > Q will be called a T-constrained inertial motion if, in any 

set of coordinates, its local expression q(-) is a solution of 

and satisfies the non holonomic constraints 

q € r. (6.3) 

Definition 6.3 Let S C TQ be any set. We say that S is invariant for the free inertial flow, 
(or equivalently: inertially invariant), if, for every free inertial motion q : X — >■ Q such that 
(q, q)(io) £ 5* for some t$ £ X, one has 

(q,q)(t) € S for all teX. (6.4) 

Definition 6.4 Let S C TQ 6e any sei. M^e say i/tai S is invariant for the T-contrained in- 
ertial flow, (or equivalently: T-inertially invariant), if, for every T -constrained inertial motion 
(see Def. 6.2) q : X — > Q such that (q, q)(to) € S for some to £ X, one has 

(q,q)(t) € S for all t EX. (6.5) 

In the next theorem we state a characterization for the T-inertially invariance of set 

{TQ) HI = rn(rnA) 1 . 

Theorem 6.1 The following conditions are equivalent: 

(1) If F = 0, for every C 1 control u(-) and every solution q(-) of the control differential 
equation (on Q) 

q = h(u), 

the map 

t i y (q(t),p(t)) = (q(t) , k(u(t)j) 
is a solution of (2.11). In particular, pi(t) = and hence qi(t) = for all t G X. 

(2) For every local chart (q) the (vector-valued) quadratic form ^ defined at (4-7) vanishes 
identically. 

(3) For every local chart (q), the (vector-valued) quadratic form v i— >■ 6[[v,v] in (4.6) van- 
ishes on the subspace Pj II (M N+M ). 

(4) The sub-bundle {TQ)jjj = V n (r PI A) 1 - is T-inertially invariant. 
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Proof. We first prove the implication (1)=>(4). Let (1) hold, and let q : 1 i-> Q be 
a T-constrained inertial motion such that qi(to) = 0. Choosing u(t) = n(q(t)), it is clear 
that the equations (6.2)-(6.3) imply (2.11), with F = 0. Hence, if (1) holds, then qi(t) = 
for all tel. By construction, the non-holonomic constraint yields qn(t) = 0. Therefore 
q(t) = qm(t) € {T q Q)jjj for all 4 € X. 

Next, we prove the implication (4)=>(2). Fix any point (7(0) = go and any vector loq € 
iR M . Let £ 1 y q(t) be the T-constrained inertial motion starting at qo, with q(0) = h(oJo). Then 
(?(•) provides a solution to the system (2.11) in connection with the control u(t) = ir(q(t)). 
If (4) holds, then qi(t) = for all t, and hence pi(t) = as well. On the other hand, our 
construction implies u(0) = ujq. Therefore, at time t = 0, the second equation in (4.8) yields 

= pj(0) = *[«(0),u(0)] = *[wo,wo]. 

Since qo and u>o can be chosen arbitrarily, we conclude that ^ = 0, proving (2). 

We now prove the implication (2)=>(1). Let a continuously differentiable control u(-) be 
given. If condition (2) holds true and F = 0, then in the ODE (4.8) for pj each term on the 
right hand side is a linear or quadratic w.r.t. pi. Therefore, if pi(to) = at some time to then 
pi(t) = for all times t G J, proving (1). 

Finally, we show that (2)<s=>(3). Since the linear map k : T n ^U i-> (T q Q)/// is a bijection, 
the equivalence of (2) and (3) follows directly by the definition of 9i and ^. □ 

Remark 6.1 In the holonomic case the the quadratic term in it, which prevents the system 
to be "fit for jumps", is completely determined by the relation between the kinetic metric 
and the distribution A. In particular this quadratic term accounts for the curvature of the 
orthogonal distribution A 1 - (see [8] for a precise result in this direction). On the other hand, 
when an additional non holonomic constraint q G T is present, the orthogonal distribution A 1 - 
can have zero curvature 2 while the dynamical equations still contain a term which is quadratic 
in ii. This is the case of the Roller Racer, our first example in Section 8. 

Given a holonomic system which is fit for jumps, Theorem 6.2 below provides a suf- 
ficient condition in order that the system remains fit for jumps also after the addition of 
non-holonomic constraints. This happens in the case of the rolling ball, our second example 
in Section 8. 

Theorem 6.2 Assume that the external force F does not depend on velocity 3 . Suppose that 
the following conditions are satisfied: 

(i) The orthogonal distribution A 1 - has zero curvature, i.e.: without the nonholonomic con- 
straint q 6 T the system would be fit for jumps. 

(ii) There exists a A-adapted system of coordinates (q) = (%<&>), with q$ = (q 1 , . . . , q N ), 
Q\> = (q N+1 , ■ ■ ■ ,q N+M ) , A = span{-£p: , . . . , ^w} , such that the Hamiltonian H is 
independent of the coordinates q\, . 

2 We recall that A x has zero curvature if the following holds. Let s q(s) be any geodesic of the manifold 
Q which intersects perpendicularly one of the leaves of the foliation generated by the integrable distribution 
A. Then every other leaf touched by this geodesic is also crossed perpendicularly [8, 24, 26, 27]. 

3 We can relax this hypothesis by allowing F to be affine in the velocity. 
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(Hi) There exists a basis {Vi, . . . ,V~n+m} forTQ such that, in this basis, the projection V\ 
is represented by a constant (i.e., independent of qj matrix Pf . 

Then, even after the addition of nonholonomic constraints q G V , the system remains fit for 
jumps . 

Proof. We need to show that the time derivative pi is an affine function of u. Recalling 
(4.5) we have 

Pi = (^■q)(p)+PI(p). (6.6) 

dP* 

By the assumption (iii), = 0. Therefore 

pi = pup), 

where p(-) solves the second equation in (2.12), namely 

dH „ „ 
p = -- + F + R. 
dq 

dH 

By the hypotheses on F and since Ri = it is then sufficient to show that — — — can be 

oq 

expressed as a function of (q,pi,u) which is affine in it. The assumption (ii) yields 

dH ( dH dH\ ( dH \ 
dq \ dq§ dq\, J \ dq t ' J ' 

On the other hand, by considering the cotangent natural basis (dq^, dq\,) and the corresponding 

dH 

adjoint coordinates (phP\,), we can express — — — as a function of {q,p*, %), 

()< b 



dH dH 



(q, Ph P\>{q^Ph^)\ 



where p\, = p\,(q,Pi, %) is the function obtained by partially inverting the linear isomorphism 

dH dH 
q = ——. By the assumption (i), the function — — — (considered as a function of the variables 
dp dq$ 

(q,Pfi,q\))) is in fact affine w.r.t. q\, (see e.g. [4] or [24]). This is equivalent to saying that it 

is affine w.r.t. it, because the map u i-> q\,(u) is a linear isomorphism 4 . Let n$ denote the 

projection p i-> n$(p) = (p$,0). By (4.3) we have 

(PB,0) = TT t (p) = ir$(pi+pni) = n$(pi + k(u)). 

Therefore — |^ is a function of (q,pi, u) which is affine in u. □ 



For a suitable choice of local coordinates on U, this isomorphism coincides with the identity. 
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7 A geometric interpretation of the quadratic term 



We begin this section by observing that it is not restrictive to assume that the manifold Q is 
an open subset of an Euclidean space. 

Indeed, by the isometric embedding theorem for Riemann manifolds [14, 21], there exists a 
positive integer M' and an embedding i : Q \— > ]R N + M + M (with ]R N + M + M endowed with the 
usual Euclidean metric) which preserves the Riemann metric on Q . Let {ei, . . . , ejv+M+M'} be 
the standard orthonormal basis of ]R N + M + M . Fix any point q €E Q. By the implicit function 
theorem, we can assume that in a neighborhood of i(q) the image of Q can be written as 

{(x 1 ,...,x N+M +M>); Xi = f i (x 1 ,...,x N+M ), N + M + l<i<N + M + M'Y 

For x = (x±, . . . , xm+n+M'), we denote by q(x) the unique point q in a neighborhood of 
q G Q such that the first N + M coordinates of i(q) coincide with (xi, . . . , xn+m)- 

Next, consider the enlarged control manifold W = U x M M . Together with the submersion 
7r : Q^-U, consider the submersion tt' : ]R N + M + M ,_>. K' by setting 

7r'(xi, . . . ,x n+ m+m>) = (u, (km+i, • • • , u M +M')) € U x 1R M ' , (7.1) 

where 

u = vr(q(x)), m = Xi- fi(x u ...,x N+M ) for N + M + 1 < i < N + M + M' . 

Calling Dl the differential of the embedding i at the point q(x), we use = Dlo A q ( x ) to 
denote the holonomic distribution on ]H N + M + M ' generated by the submersion (7.1). Given 

the non-holonomic distribution V on Q, the corresponding non-holonomic distribution V on 
M N +M +M> ig defined ag 

T x = span^DL o r q(x) , ejv+M+i, ••• ,e N+M +M'}- (7.2) 
Observe that, if (2.6) holds, then we also have 

A x + r x = m n+m+m '. 

Trajectories of the non-holonomic system with active constraints defined at (2.1)-(2.2) can now 
be recovered as trajectories of the non-holonomic system with active constraints (7.1)-(7.2) 
corresponding to control functions 

t ^ (u(t),u M +i{t),...,u M +M'(t)) e UxM M ' 

with UM+i(t) = ■ ■ ■ = UM+M'(t) = for every time t. 

Based on the previous remarks, without loss of generality we now analyze the geometric 
meaning of the term ^[u, it] in (4.7), assuming that Q = M N+M is a finite dimensional 
vector space with Euclidean metric. In this case the coefficients of the metric are constant: 
grs = 9 rs = 5rs- Moreover, with canonical identifications we have 

P = q, Hu) = h(u), %- ^ Pi = Pi- (7-3) 

Therefore, (4.6)-(4.7) yield the simpler formula 

*[«,*] = (^■Ku)){Ku)). (7.4) 
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In this setting, consider a system of local coordinates (u±, . . . , um) on U. To fix the ideas, 
let v = (1, . . . , 0) € M M . We seek a geometric characterization of the term 

nv,v] = • %))(/»(«)) • (7.5) 

At each point g e JR M+7V in a neighborhood of choose an orthogonal basis of unit vectors 
Vi(g), . . . , Vtv+m('z) as in (5.1), but with V^+i being the unit vector parallel to h(v), so that 

h(v,q) = £(q)V N+ i(q) 

for some £(q) > 0. Notice that we now write h = h(v,q), since we regard h as a vector field 
on Q = M N+M . In the following, we use (•, •} to denote the Euclidean inner product, and we 
write (D(j))'Vi for the directional derivative of 4> in the direction of Vj. Moreover, the map 
t i y (exptVj)(x) denotes the solution of the Cauchy problem 

x(t) = Vi(x(t)), x(0) = x. 

Computing (7.5) in terms of the basis {Vi, . . . , Vn + m} we obtain 

N-u 

? N+M 



Pj(w) = J2\ w ' V V V * fora11 w^M r 
i=i 

N-u N-u 

«] = E (^JV+i , (£>Vi) ?V W+1 ) V, + E (eVjv+i , (DV^Vat+i 



A-l t=1 (7-6) 

= E (CViv+i, (^V^Vjv+i)^. 
1=1 

We now claim that, for 1 < i < N — v, the Lie bracket 

[ey N+l , V,] = (DV i )ZV N+1 - (D(£Vjv+i))Vi 
lies in the subspace A q and hence is orthogonal to Vjv+i- Indeed, for every q one has 

[ev w+ i, vm = iim ^t^, 

where 

q(e) = (exp(-eVi)) o (exp(-e£VAr + i)) o (exp(eVj)) o (exp(e£V;v+i)) (q) . 
However, for any q we have 

-^(expftVi))^) = 0, ^7r(exp(teV w+1 ))(g) = u . 

Hence, for every e one has 

7r (<?( £ )) = (exp(-£v)) o (ex.p(£vf)(Tr(q)) = n(q) . 

This implies 

[€Vn+i, Vi](g) G kerDir(q) = A q , 
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as claimed. Therefore 

(v N+1 ,(DV i )(€V N+1 )) = (Vjv+i, {D(€V N+1 ))Vi) . (7.7) 
Next, we observe that 

(Vat+i, (^Vjv+OVi) = ±(D\V N+1 \ 2 )Vi = 0, 
because all vectors Vjv+i(g) have unit length. From (7.6) and (7.7) we obtain 

N—v 

*[v,v] = (CViv+i, (£>(^ + i))Vi)Vi 

i=i 

(7.8) 

i=l 

Recalling (3.14), we can think of £(<?) = \h(v,q)\ as the infinitesimal distance between 
leaves of the foliation generated by the submersion it at the point q, computed along paths 
that respect the non-holonomic distribution. According to (7.8), if the directional derivative 
of this distance is zero in all directions of the subspace (T q Q)j = A q nT q , then the term 
^[v, v] vanishes for all v G T n ^U. 

Going back to the general case of a Riemann manifold Q and a submersion it : Q i— >■ U, we 
thus obtain the following geometric characterization. For every u G U, consider the leaf 

A u = {qeQ; vr(q) = u}. (7.9) 

Then for every q G A u and every tangent vector v G T U U there exists a unique vector 
h(v) G (TqQ)jjj such that Dtt ■ h(v) = v. Recalling (3.14), one could define an inner product 
on (•, -} u on T U U by setting 

(v,v) u = g q [h(v), h(v)]. (7.10) 

Of course, in general this inner product depends on the choice of the point q along the leaf 
A u . According to (7.8), the condition ^ = means that the directional derivative of the right 
hand side of (7.10) vanishes along A q n T q . More precisely: 

Proposition 7.1 The following conditions are equivalent. 

(i) * = 0. 

(ii) For every u G U, v G T U U, and every smooth path s h-> q(s) G A u such that q(s) G 
A q ( s ) nr q ( s ), the inner product g q ( s )[h(v), h(v)] is constant. 

Remark 7.1 In general, the condition ^ = does not guarantee the existence of a Ricman- 
nian metric on U defined by (7.10). In order that this metric be well defined, i.e. independent 
of the point q G A u , one also needs to assume that the distribution A q n T q is completely 
non-integrable restricted to each leaf A u . More precisely, every two points qi,q2 G A u should 
be connected by a path t H> q(t) G A u with q(t) G A q( - t ) n r q ( t ) for all t. Notice that this 
condition is clearly satisfied if the non-holonomic constraint is not present (i.e., T q = T q Q) 
and the leaves A u are connected. From (7.8) we thus recover some earlier characterizations of 
the property "fit for jumps" given in [9, 24]. 
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Figure 2: The roller racer. 



8 Two examples 
8.1 The Roller Racer 

As a toy model, consider the Roller Racer, which will serve as a simple illustration of the 
general theory. This is a classical example of a non-holonomic system, widely investigated 
within the theory of the momentum map [2, 12]. It consists of two rigid planar bodies, 
connected at a point by a rotating joint, as shown in Figure 2. One of the two bodies has a 
much larger mass than the other. Let p be the distance between the joint and the center of 
mass of the heavier body. To simplify computations we also assume that the center of mass 
of the lighter body coincides with the joint. 

The coordinates (q 1 ,q 2 ,q 3 1 u) are as follows . We let (<j ,3 ,(7 1 ) = (x,y) be the Euclidean 
coordinates of the center of mass of the large body. Moreover, q 2 = 9 is the counter-clockwise 
angle between the horizontal axis and the major axis of the large body, and u = 4> is the 
counter-clockwise angle between the major axes of the two bodies. 

The non-holonomic constraint consists in assuming that each pair of wheels roll without 
slipping, parallel to the corresponding body. This corresponds to the condition 



In this case the (non integrable) distribution is 2-dimensional. The admissible motions t i— > 
(q,u)(t) are thus subject to 



(q,u) G T = ker(oj 1 )Dker(uj 2 ), 



where 




(8.1) 



cosq 2 ^)^^) -sin q 2 (t)q 3 (t) = 0, 



< 



(8.2) 



cos(q 2 (t) +u(t))q 1 (t) +p cos u(t) q 2 (t) - sin(g 2 (t) +u(t))q 3 (t) = 0. 
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The coordinate (f) = u is regarded as a control. Notice that the transversality condition 
(2.6) is trivially satisfied, because 

A ker = span{du}, T ker = spaniuj 1 ,lo 2 }, 

and spanluo 1 , lo 2 } n span{du} = {0}. 

For simplicity, we assume that the large body has unit mass. The moment of inertia of 
the large body w.r.t. the vertical axis passing through its center of mass is denoted by I. The 
mass of the small body is regarded as negligible, but its moment of inertia J w.r.t. the vertical 
axis passing through the center of mass (of the small body) is assumed to be different from 
zero. 5 The kinetic matrix G = (<?ij) and its inverse G _1 = (g l 'i) are computed as 



(9, 
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Let us start by finding a basis for the decomposition TQ = {TQ)j © (TQ) n (TQ) in . 
Notice that the vectors 



o d . . d 2 d 

Wi = 2p cos u sin q -— + 2 sin it— ^ + zpcosucosq 
oq l oq z 



dq 3 ' 



w 2 



d_ 

du 



form a basis for T. Since A = span{-^, -J^, -£js}, {wj is a basis for (Tq)j = T n A. 
It is also straightforward to verify that 

(TQ) n = T x = ker{g( Wl ),g(w 2 )} 

where 

/ esc q 2 tai 

V 2 = 7TT — TT^T + ~7T~ j V 3 



dq 1 dq 2 du' 



span{v 2 ,v 3 } , 

cot 2 — ^ — 

^ dq 1 dq 3 



(8.3) 



Since (TQ) in is orthogonal to (TQ)/ © {TQ) n , 

{TQ)m = ker |g(wi),g(v 2 ),g(v 3 )| = span{v 4 } 



(8.4) 



with 



v 4 



1 

A^ 



1 _ , 2 • r • 2 ^ 1 t 2 ■ ^ 3 

-Jpsmq sm2ii — — r — J sm 1x77-77 Jpcosq sm2n — 

2 ag 1 ag^ 2 dq 6 



+ 



d_ 

du ' 



we have Ao = p 2 cos 2 u + (I + J) sin 2 it. 

Defining Vi = wi,V 2 = v 2 , V3 = V3, V4 = V4 we obtain a family of vector fields (on an 
open subset of the configuration manifold) as in (5.1). 

To obtain the equations of motion we need to compute the coefficients (5.9). For this 
purpose, let us begin with observing that the form fl\ generating (T*Q)j is given by 

Ol = n^dq 1 + n 2 ,ldq 2 + n 3A dq 3 + n 4A dq 4 = g(Vi) 

= (2p cos u sin q 2 )dq 1 + (2(1 + J) sin u)dq 2 + (2p cos u cos q 2 )dq 3 + (2 J sin u)dq A . 



'This is a standard approximation adopted in the existing literature. 
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Let us recall that the projection matrix Pf is computed as (see (5.7)) 



(P*Y = V g r ' 

fc=l 2^a,b=l9a,bV I* 1 

Let us set Ai = / + J + p 2 + (—7 — J + p 2 ) cos 2u. Then an elementary computation of (5.8) 
yields the following system of four differential equations: 



q 1 = 


2p cos using 2 • £ 


Jp sing 2 sin 2u 

2A U ' 


e = 


2sinu-e Js i n2 




q 3 = 


2pcos q 2 cos-u ■ £ 


Jp cos g 2 sin 2u „•, 
2A U ' 


£ = 


2(7+ J-p 2 ) sin 2m 
Ai 


■£u+ 8J P 2 A ? su -u 2 



(8.5) 



Remark 8.1 Assume that we implement a vibrating control, of the form u(t) = u+eK sm(t/e) 
for K G K some e > small. By the results from [7, 8] , each solution of the system 

2pcosu smq 2 ■ £ 
2 sin u ■ £ 
2p cosu cos q 2 ■ £ 
4 J/j 2 cos u • 7C 2 



(8.6) 



A? 



(formally obtained by neglecting higher order terms and then averaging), can be uniformly 
approximated by the solutions of the original system 8.5 6 . It is then easy to verify a kine- 
matical well-known behavior of the Roller Racer: fast oscillations of the handlebar around a 
angle u produce a motion which, asymptotically, is a superposition of forward motion and a 
rotation of the angle q 2 (which is the angle between the principal axis of the large body and 
the x-axis, namely the q 3 axis). In particular if the oscillations are around the central position 
u = 0, then q 2 is constant and the motion approaches asymptotically a forward motion. Of 
course this agrees with results obtained with methods based on momentum maps, where, in 
particular, controls are forces or torques. See [2] or the extensive analysis of the Roller Racer 
in [15]. 



8.2 A ball rolling on rotating disc 

The next example illustrates an application of Theorem 6.2. It represents a modification of 
the classical case (see [2] ) of a mechanical system consisting of a ball that rolls without sliding 
on a disc which rotates with constant speed (Figure 3). The somehow surprising fact is that 
the ball does not move away toward infinity, regardless of the angular speed of the disc. In 

6 More generally, one could approximate any solution of a (possibly) impulsive system obtained 8.5 by 
replacing (ii,u 2 ) with a pair (a, fx), where a is a L 1 control and /i is a positive measure such that a 2 < n (in 
the sense of measures) 
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fact, the barycenter of the ball remains confined within a bounded region, depending on the 
initial conditions. The "surprise" actually comes from the fact that intuition suggests a sort 
of centrifugal effect, pushing the ball outward in the radial direction. 

We slightly modify this example by allowing the rotational speed of the disc to be not 
constant. More precisely, we consider the rotational angle u of the disc as a (state-space) 
parameter u, regarded as a control. We will show that, as in the case of constant rotational 
speed, a small rapid oscillation of the disc around the center does not instantly push the ball 
toward infinity. As a consequence, the system is "fit for jumps" . To recast this system in our 
general framework, consider a material disc D, centered at a point O with inertial moment 
/ > 0. If u(t) designates the counter-clockwise angle of the disc with respect to inertial 
coordinate axes, its kinetic energy is given by \lu 2 ■ Actually, in our setting the value of I 
is irrelevant for the motion of the system, because the angle u is taken as a control variable. 
The configuration manifold of the rotating disc together with the rolling ball can be identified 
with Q = M 2 x SO(3) x S 1 . The first two coordinates (x, y) € M 2 describe the position of the 
contact point w.r.t the inertial frame. A unitary matrix R € 50(3) represents the rotation 
of the ball, while the angle u G S 1 describes the rotation of the disc. Here the submersion 
describing the control constraint is simply the projection 

7r : Q — > S 1 , ir(x, y, R, u) = u. 

In addition, the non-holonomic constraints are represented by the linear equations 

x + roj 2 + uy = , 

(8.7) 

y — ru\ — ux = , 

where r is the radius of the ball and oo = (001,002,003) is the angular velocity vector of the 
ball w.r.t. the inertial frame. For simplicity, we assume that the inertial moment of the disc 
w.r.t. its center is I = 1 and that the ball has homogeneous density and total mass equal to 
1. The kinetic energy T of the system is then given by 

T = ijW^ + ^E^ + u 2 ), (8.8) 

where k 2 is the moment of inertia of the ball w.r.t. any of its axes. 

Let (q 1 , . . . , (/ 6 ) be local coordinates for Q, where q 4 = x, q 5 = y, c/ 6 = u, while (q 1 , q 2 , q s ) 
are local coordinates for SO(3). Recalling the meaning of angular velocity 00 and the form 
of the kinetic energy T, we can find an orthonormal basis {Vi, . . . , Vg} of the tangent space 
TQ such that the following holds: 

d 

Vi = g-r for ie {4,5,6}, 

v, = E4(^V)^ for {1,2,3}. (8.9) 

Moreover, if to = (001,002,003) describes the angular velocity of the ball, then 

3 

(91,92,93) = E KW i v j- 
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Figure 3: A ball rolling without sliding on a rotating disc. 



In terms of these coordinates, if q = J2j c j^ji the non-holonomic constraints (8.7) take the 
form 

Q + -C2 + q 5 c 6 = , c 5 c\ - q A c 6 = . 

K K 

We thus obtain 

A = span{Vi,V2,V 3 ,V4,V 5 }, A ± = span{V 6 }, (8.10) 

ART = span (vi + -V 5 , V 2 - -V 4 , V 3 ) , (8.11) 
L re re J 

(Afir) 1 = spanjVi-^Vs, V 2 + ^V 4 , V 6 J. (8.12) 

Using the basis{Vi, . . . , Vg}, all projections Pj,Pj, J € {1,11,111} are given by constant 
matrices, so hypothesis (hi) in Theorem 6.2 is verified. By (8.8)-(8.10) it is clear that the 
matrix g of the kinetic energy with respect to the basis (^t, • • • , -^) is independent of q e and 
has all zeros in the sixth row and the six column, except for the corner entry which is equal 
to 1. Therefore the same is valid for the inverse matrix g~ l . Denoting by p T the transpose of 
the row vector p, this implies 

• H = \pg~ x p T is independent of q 6 . 

• If the non-holonomic constraint (8.7) is not present, then the system is fit for jumps. 

By Theorem 6.2 we thus conclude that this system is fit for jumps. In particular, a rapid 
small oscillation of the disc around its center will not cause the ball to fly away to infinity. 
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